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Graphene - a two-dimensional (2D) material with 
unique electronic properties appears to be an ideal object 
for the application of surface science methods. Among 
them, a family of scanning probe microscopy meth¬ 
ods (STM, AFM, KPFM) and the corresponding spec¬ 
troscopy add-ons provide information about structure 
and electronic properties of graphene on the local scale 
(from /rm to atoms). This review focuses on the recent 
applications of these microscopic/spectroscopic meth¬ 
ods for the investigation of graphene on metals (inter¬ 
faces, intercalation-like systems, graphene nanoribbons 
and quantum dots, etc). It is shown that very important 
information about interaction strength at the graphene- 
metal interfaces as well as about modification of the elec¬ 
tronic spectrum of graphene at the Fermi level can be 
obtained on the local scale. The combination of these re¬ 
sults with those obtained by other methods and compari¬ 
son with recent theoretical data demonstrate the power of 
this approach for the investigation of the graphene-based 
systems. 



STM image of a graphene island (quantum dot) on Ir(l 11). 
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1 Introduction The discovery of the fascinating elec¬ 
tronic transport properties of graphene EIIIEI, a two- 
dimensional (2D) allotrope form of carbon, has attracted 
intense attention of materials science and solid state physi¬ 
cists ii4im . because of fundamental physics as well as 
the device application perspectives in electronics and 
spintronics. First exfoliated and identified in 2004 |fni2l 
a, graphene was actually known for many years in the 
surface science community as a graphitic “dead” layer, 
which poisons the catalitic activity of metal surfaces El 


m[8ll9l . However, the interest in the graphene/metal sys¬ 
tems experienced a strong revival since it has been realized 
that the synthesis of graphene on metal surfaces is the 
most prospective method to obtain high quality large-area 
graphene samples for further applications nniiniini. 

Beyond the possible technological relevance, the inves¬ 
tigation of continuous graphene layers as well as graphene 
nanostructures on metal surfaces turned out to be chal¬ 
lenging and exciting from both experimental and theo¬ 
retical point of view. One of the most prominent exam- 
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Figure 1 (a) Principal scheme of an SPM experiment, (b) 
Energy scheme of the STM experiment for different bias 
voltages applied between the tip and the sample, (c) Fre¬ 
quency shift of the sensor (Af) as a function of the dis¬ 
tance between tip and sample (z) expressed from the to¬ 
tal interaction energy (Etot)', A: no interaction, B; long- 
range vdW and electrostatic forces, C: short-range chem¬ 
ical forces giving atomic contrast in AFM, D: short-range 
repulsive interaction. 


pies in this regard addresses the nature of bonding at the 
graphene/metal interface, the problem, which was dis¬ 
cussed in a large series of publications (see e. g. Refs. uni 
[TrlfTSl l. having the question open for long time - how rel¬ 
atively weak interaction (in the order of 50 — 150meV/C- 
atom) can lead to strong modihcations in the electronic 
structure of graphene? Only recently this problem was 
considered in a complex approach and the universal model, 
describing the interaction in the whole graphene/metal 
interface family, was proposed M- The search for the 
routes to minimize this interaction especially aims at the 
preparation of the graphene nanoribbons 017lll8lfT9llMl 
l2ni22l or epitaxial nano sized islands with different edge 
terminations ||23|24l|25l|26l, where the low-dimensional 
effects such as a band-gap width quantization and the 
edge-induced magnetism are expected. Due to the truly 
2D character of the crystallographic and the electronic 
structure of graphene and the localisation of the inter¬ 
esting phenomena at the small scale, the scanning probe 
microscopy (SPM) methods, i.e. scanning tunnelling mi¬ 
croscopy (STM) and atomic force microscopy (AFM) in 
combination with the corresponding spectroscopy tech¬ 
niques (STS and AFS), can provide important information 
about the properties of graphene nanostructures at the 
nanometer- and atomic-scale. 

This feature article gives a short review of the recent 
STM and AFM experiments on continuous graphene layers 
and graphene nanostructures prepared on metallic surfaces. 
The presented results are linked to the available theoretical 
data providing an important input for the understanding of 
the observed phenomena. 

2 Scanning probe microscopy and spectroscopy 


2.1 STIUI and STS In STM experiments the sharp 
conductive tip is moved above the sample surface in x 
and y directions [Fig. [TJa)]. In constant current STM 
measurements the tunnelling current It, which exponen¬ 
tially depends on the distance z between tip and sam¬ 
ple [It oc exp{—kz)], is kept constant by a feedback 
loop adjusting the actual z position, which is recorded 
to obtain a three-dimensional topography map. The ob¬ 
tained constant current image represents the integral 
local density of states (FDOS) of the sample surface, 
It{x, y, z,V) oc f p{x, y, z, E)dE, where the integration 
is performed in the energy range between the Fermi lev¬ 
els {Ep) of tip and sample, when they differ by the value 
cUt, where Ut is the bias voltage [Fig. [TJb)]. Differenti¬ 
ation of the / [x, y, z, Ut) with respect to the bias voltage 
(can be performed with lock-in technique) gives a direct 
information about the local density of states. Such mea¬ 
surements carried out in the scanning mode at different 
bias voltages allow to observe the so-called electron den¬ 
sity standing waves of different periodicities that can be 
used to obtained the electron dispersion relation, E{k), of 
the surface electronic states. 

2.2 AFM and AFS In the modern AFM experiments, 
frequency-modulated AFM (FM-AFM), the oscillating 
at high resonance frequency scanning sensors (based on 
quartz tuning forks or length-extensional resonators) are 
widely used II27II28L In this case the sensitivity of the 
method can be increased dramatically. Here, a tiny con¬ 
ductive tip connected to the resonator approaches the 
sample surface and the interaction between them leads 
to the change of the resonance frequency by value Af: 
F^{d) = dE/dd, Af{d) = -/o/2fco • dF^{d)/dd, where 
E{d) and Fz{d) are the interaction energy and the vertical 
force between a tip and the sample, respectively, /o and 
fco are the resonance frequency and the spring constant of 
the sensor [Fig.[T][c)]. The detected frequency shift, Af, is 
used in the feedback loop and the corresponding “topog¬ 
raphy” map can be obtained in such AFM measurements. 
The resulting interaction between the sample and the tip is 
the sum of the long-range electrostatic (Fei) and van der 
Waals {Fydw) contributions and the short-range chemi¬ 
cal (Fchem) interaction. Similar to the STM experiments, 
Fchem at short distances between tip and sample produces 
the atomic contrast in AFM experiments and as was shown 
in a series of the recent works the observed imaging con¬ 
trast depends on several factors, e.g., distance and the 
setting point for the measurements [Fig. [TJc)]. The locally 
measured F^ (d) curves can be used for the calculation of 
the tip-sample interaction forces via methods developed 
by Giessibl li^ or Sader and Jarvis OOl . 

As an add-on to AFM, the method of the Kelvin-probe 
force microscopy (KPFM) allowing to obtain the local dis¬ 
tribution of the electrostatic potential was developed m 
[^ . In this method the topographic measurements are per¬ 
formed in AFM mode and at the same time the DC (Udc) 
and AC voltages are applied between tip and sample. Here 
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(a) Top View (c)(1»1) structures 



Figure 2 (a,b) Top and side views of the crystallographic 
model of the graphene moire structure on the close- 
packed (111) metallic surface, (10 x 10)graphene/(9 x 9) 
Metal(l 11). (c) Structures of the local high-symmetry po¬ 
sitions of the graphene/metal interfaces. 

the lock-in technique is used and the hrst harmonic in 
the output signal, which is proportional to the difference 
between Udc the local contact potential difference 
(Ulcpd) between conductive tip and sample, is nullihed 
allowing to extract the ULCPD{x,y) map of the sample 
surface. 

2.3 3D AFM/STM The 3D STM/AFM measurements 
can be performed in two ways. In the hrst case the dis¬ 
tance between tip and sample is hxed with the feedback 
loop completely switched off and the respective I{x,y) 
and Af{x, y) maps are simultaneously collected. If the 
distance between tip and sample is varied with the reg¬ 
ular steps, then the dense 3D data sets can be collected, 

I{x, y, z) and Af{x, y, z). In the second approach the sin¬ 
gle I(z) and Af{z) curves are measured on the {x, j/)-grid 
and then they are combined in the 3D sets. Both measure¬ 
ment schemes are time consuming and here the special at¬ 
tention should be paid to the possible creep and thermal 
drifts of the piezo-drive and the AFM sensor as well as to 
the post-correction of the obtained data. Different schemes 
of 3D STM/AFM measurements and their (dis)advantages 
are discussed in Refs. 033ll34l[T5l[3^ . 

3 Theoretical approaches: DFT, models, etc. 

Considering the arrangement of the graphene layer on a 
metal surface, whether lattice-matched or lattice-mismatched 
system, one can identify several high-symmetry stacking 
positions for carbon atoms in the layer. (Note: Different 
notation are used in the literature to mark these positions). 
They are: 

(a) ATOP (hcp-fcc) position, where carbon atoms surround 

the metal atom of the top layer and are placed in the hep 

and/cc hollow sites of the Metal(lll) stack above (S- 

1) and (S-2) Metal-layers, respectively [Circle in Fig.[^ 


(a) and the respective arrangements in Fig. (b) and 

(c)]; 

(b) FCC {top-hep) position, where carbon atoms surround 
the/cc hollow site of the Metal(lll) surface and are 
placed in the top and hep hollow positions of the 
Metal(lll) stack above (S) and (S-1) Metal-layers, 
respectively [Rhombus in Fig.|^(a) and the respective 
arrangements in Fig.[^(b) and (c)]; 

(c) HCP {top-fee) position, where carbon atoms surround 
the hep hollow site of the Metal(lll) surface and 
are placed in the top and fee hollow positions of the 
Metal(lll) stack above (S) and (S-2) Metal-layers, 
respectively [Triangle in Fig. (a) and the respective 
arrangements in Fig.|^(b) and (c)]; 

(d) BRIDGE-position, where carbon atoms are bridged by 
the Metal atom of the (S) layer [Rectangle in Fig.|^(a) 
and the respective arrangements in Fig.|^(b) and (c)]. 

The vast majority of computational studies of graphene 
on metals are currently performed using hrst-principles 
electronic structure methods based on density functional 
theory (DFT). While density functional theory itself is ca¬ 
pable of providing the exact solution to the Schrodinger 
equation, including long-range correlations - the disper¬ 
sion, the approximations made in the DFT functionals 
of all types yield a rather unsatisfactory description of 
the intermolecular interactions. Therefore when studying 
graphene/metal interfaces, van-der-Waals bound systems, 
one usually employs a posteriori dispersion correction 
schemes, such as DFT-D method of Grimme 037113811^ or 
Tkachenko-Scheffler method Il40l . and the obtained results 
are usually in reasonable agreement with experiment. Un¬ 
like these post-DFT corrections, an a priori way consists in 
utilisation of non-local correlation functionals that approx¬ 
imately accounts for dispersion interactions, the so-called 
vdW-DF Il4ni42l . Although the latter approach is expected 
to give more accurate results, it is shown to predict irra¬ 
tionally weak binding of graphene on metals surfaces for 
all studied cases li43l . 

A very important analysis of the surface atomic and 
electronic structure is the simulation of STM images. 
Among early theoretical models to describe atomically 
resolved STM images Tersoff and Hamann approach ll44ll 
has found particular attention since it relates the observed 
contrast to a conceptually simple quantity of the surface, 
namely the local density of states of valence or conduc¬ 
tion electrons. While it provides a simple interpretation 
of the STM results, this model ignores the influence of 
the tunnelling barrier, the structural and electronic proper¬ 
ties of the tip and its interaction with the sample surface. 
However, more exact theoretical formulations of the STM 
problem (for a review, see Ref. USD were found to be 
much more time-consuming. At the same time, in all stud¬ 
ied graphene/metal cases the Tersoff-Hamann approach 
was found to be rather reliable. 

While STM can achieve atomic resolution, it probes the 
surface electronic structure of the sample in a mode that 
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Figure 3 STM images of (a) lattice matched 
graphene/Ni(lll) (2.5 x 2.5nm^, Ut — 2mV, 
It = 48nA) and lattice mismatched systems; (b) 

graphene/Rh(lll) (5.4 x 5.4nm^, Ut = 300mV, 

It = 1.6nA), (c) graphene/Ru(0001) (5.7 x 5.7nm^, 
Ut = 300mV, It = 1.6nA), and (d) graphene/Ir( 111) 
(4.5 X 4.5nm2, Ut = 300mV, It = 1.6nA). 


does not provide a direct interpretation of the atomic struc¬ 
ture. First-principles simulations of NC-AFM can enhance 
the interpretation of experimental measurements; however, 
such simulations remain a challenge because they involve 
calculations of the sample together with an atomic model 
of the AFM tip (see e.g. Ref Il46l ). Hence, the simulation 
can be very laborious. An efficient scheme to simulate NC- 
AFM images using inputs from first-principles calculations 
of the sample only without explicit modelling of the AFM 
tip was proposed by Chan et ai ll47l . Both of these ap¬ 
proaches were found to work reasonably well for the con¬ 
sidered graphene/metal systems. 

4 Graphene on metals: SPM/DFT of lattice match¬ 
ed and lattice-mismatched systems The history of 
the surface science studies of graphene on metals can be 
traced back to the middle of the 60s. However, the first 
STM experiments on graphitic layers on metal surfaces 
were performed in the beginning of the 90s when graphitic 
layers on Ni(lll) ll48l and Pt(lll) ||9l were studied. It is 
worth to mention that already at that time the graphene 
moire structures of different periodicities were identified 
onPt(lll). 

The simplest, from the structural point of view, 
graphene-metal interface is graphene/Ni(l 11) or graphene/ 
Co(OOOl) as the difference between lattice constants of 
these surfaces and graphene is below 2%. Both systems 
were intensively studied by means of STM 1131501151] 
l52ll53ll54ll55ll56lB7ll58ll . These studies demonstrate that a 
graphene layer deposited by means of chemical vapour 
deposition (CVD) from hydrocarbons at optimal experi¬ 


mental conditions forms a commensurate (1x1) structure 
on top of Ni and Co close-packed surfaces [Fig. [^a)]. In 
this case, according to DFT calculations (with and without 
inclusion of the dispersive vdW interaction) 05911 60 IIMI 
[54l . graphene is adsorbed in the HCP (top-fcc) configura¬ 
tion on Ni(lll) or Co(OOOl). However, as shown within 
the same theoretical calculations, the energy difference be¬ 
tween all high-symmetry stackings is rather small (below 
100 meV) that can lead to the appearance of the moire-like 
structures even for these lattice-matched interfaces 15^ 
|57l. 


The interesting class of graphene-metal systems is ob¬ 
tained when graphene is prepared on the close-packed sur¬ 
faces of 4fi or 5d metals. In this case, the existing lattice 
mismatch leads to the formation of the so-called moire 
graphene-metal structures of different periodicities. These 
systems were in the focus of the intensive STM inves¬ 
tigations and the representative examples are shown in 
Fig. I^b-d). It was found that the apparent corrugation of 
graphene-moire strongly depends on the imaging bias volt¬ 
age 062ll63ll64ll65ll66llff7l and at some tunnelling condi¬ 
tions inversion of the imaging contrast is observed. 

The representative example is shown in Fig. Qa) 
where bias voltage was changed during STM imaging 
of graphene/Ir(lll) from Ut = —0.3V (bottom) to 
Ut = —1.8V (top) 067II68I . Here at the typical tunnelling 
voltages used in STM, the graphene/Ir(lll) system is 
imaged in the so-called inverted contrast when topograph¬ 
ically highest places of the moire structure (ATOP=A) are 
imaged as dark areas and the lowest places are images as 
bight (FCC=F and HCP=H) [Fig. I^^a), bottom]. Beyond 
the change of the imaging contrast, the apparent corru¬ 
gation is also changed from « 0.32 A (Ut = —0.3 V) 
to « 0.14 A (Ut = —1.8 V). This effect was success¬ 
fully explained theoretically where graphene/Ir(l 11) was 
treated on the DFT-D2 level. These calculations correctly 
describe the crystallographic structure of graphene/Ir(l 11) 
with the corrugation of graphene of 31pm (PBE-D2) IIFtII 
or 35 pm (vdW-DFT) |[69l and reproduce the dependence 
of the imaging contrast as well as the corrugation in STM 
experiment on the bias voltage l67l . The effect of the 
contrast inversion was assigned to the so-called interface 
states of graphene/Ir( 111) formed as a result of the over¬ 
lap of graphene and Ir orbitals. The value of 
the graphene corrugation of 43 ± 9 pm was obtained in the 
LEED-IV experiments supported by the AEM results (with 
CO-terminated tip), which yield a value of corrugation of 
42 — 56pm llTOl . Contrary to these results, the much larger 
corrugation of O.GA and I.OA for 0.39ML and 0.63ML 
graphene on Ir(l 11), respectively, was obtained from an x- 
ray standing wave (XSW) experiments li^ and this large 
corrugation was assigned to the increasing of the number 
of lines of wrinkles on the surface. The similar approach 
was used for the explanation of the effect observed during 
imaging of graphene/Ru(0001) 1163 11641|65]| . 
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Lotujfit profile (nm) 


Figure 4 (a,b) Experimental (6 x 6nm^) and theoretical 
STM images for graphene/Ir(lll), respectively. In the ex¬ 
perimental data the bias voltage was changed from —1.8V 
(top) to —0.3 V (bottom) during scanning. In the simulated 
images the respective integration energy ranges between 
Ep and cUt were used, (c) Combined STM/AFM imag¬ 
ing of graphene/Ir( 111) when scanning mode was changed 
“on-the-fly” during scanning (5 x 5nm^, Up — 30mV, 
It = InA, Af = —475mHz). The lower panel shows the 
corresponding height prohles for STM and AFM. 


One of the most interesting examples of graphene 
growth on a id metal surface, which also delivers an un- 




Relatlve distance, d (pm) 



Af(x,y,d =50 pm) I(x,y.d= 50 pm) 


Figure 5 Upper panel; Af{z) and I{z) curves measured 
for different high-symmetry positions of graphene/Ir( 111). 
Lower panel: simultaneously collected pure constant 
height Af{x,y) and I{x,y) maps (3.9 x 3.9nm^) at the 
relative distance between tip and sample of d = 50 pm. 
The bias voltage was changed from -1-50 mV to —50 mV in 
the middle of the scanning area. 


usual moire structure symmetry, is the graphene/Fe(l 10) 
system GD- The preparation of graphene on top of a 
Fe(llO) him, although being challenging due to the pos¬ 
sible formation of different carbidic phases, has been suc¬ 
cessfully realized by a standard preparation method relying 
on thermal dissociation of the CaHg precursor at the hot 
metal surface. Corresponding STM images of the graphene 
surface show a pronounced corrugation yielding a periodic 
wavy pattern, with the lateral distance of 4 nm between the 
waves running along the [001] direction of the underlying 
substrate. This very special moire structure stems from a 
unique combination of the lattice mismatch at the interface 
and the strong graphene/substrate interaction as shown by 
the accompanying DFT calculations Gil. 

Application of the fast and reliable scanning probe sen¬ 
sors (Kolibrisensor™from SPECS or QPlus from Omi- 
cron) allows to perform simultaneous STM/AFM imaging 
of the same sample area and delivers interesting results 
for the graphene-metal moire structures. In AFM the in¬ 
teraction between tip and sample is the sum of the long- 
range vdW and electrostatic forces and the short-range 
chemical forces. The latter allows to perform the atomi¬ 
cally resolved imaging in AFM. In this case the imaging 
contrast in AFM depends on the distance between tip and 
sample and the corresponding Af value used in the feed¬ 
back loop of SPM. In the recent experiments performed 
on graphene/Ir(l 11) the imaging mode was changed be¬ 
tween STM and AFM (in the attractive regime) during 
scanning that allows to carefully trace the topographic con¬ 
trast in this system [Fig.|^c)]. One can clearly see that the 
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imaging contrast between CC STM and CFS AFM is in¬ 
verted; STM demonstrates inverted contrast and AFM per¬ 
formed in the attractive regime far from the minimum point 
for the Af signal gives true direct contrast. The corre¬ 
sponding corrugation of the moire structure is 100 pm and 
24 pm in STM and AFM, respectively. At the same time 
the constant-height measurements (Fig. performed for 
the repulsive part of the interaction curve of the tip and 
sample {d — 50 pm) show the stronger repulsive interac¬ 
tion between ATOP place of graphene/Ir( 111) and the tip 
compared to other regions of this structure. Change of the 
sign of the bias voltage during such measurements leads 
to the inversion of the sign of the tunnelling current (value 
of the current depends on the local density of states in the 
particular energy range around Ep), but it does not influ¬ 
ence the Af channel as the electrostatic force between tip 
and sample depends as F^i oc — C/|. and no “cross-talk” 
between Af and I channels is detected. 

The 3D force and current spectroscopy, Af{z) and 
I{z), performed on the grid for graphene/lr(l 11) demon¬ 
strates the inversion of the imaging contrast for the Af 
channel in the constant height AFM measurements m 
[68117^ . The results of such measurements are presented in 
Fig.|^a) where on the left-hand side the Af{z) curves ac¬ 
quired at the different places of the graphene moire struc¬ 
ture on Ir(lll) are shown. One can clearly see that sev¬ 
eral so-called crossing points for the Af curves appear on 
this plot that indicates that inversion of the imaging con¬ 
trast in CH AFM will be observed on the moire-cell scale 
as well as on the atomic scale. The corresponding cuts, 
Af{x, y, d), extracted from the 3D data sets are shown on 
the left-hand side of Fig. a), where such effect of the con¬ 
trast inversion as a function of the distance between oscil¬ 
lating tip and sample is clearly observed II67II68II . 

The interesting effect was observed in the AFM imag¬ 
ing and spectroscopy of graphene/lr(l 11) when clean 
metallic and a CO-passivated tip were used in the experi¬ 
ment [Fig.|^b,c)] 11721 . Imaging with the reactive metallic 
Ir tip shows that atomic contrast on graphene/Ir(l 11) can 
be obtained in the attractive regime, which is changed 
to the repulsive atomic imaging at the shorter distances. 
Here the inversion of the imaging contrast is observed 
at the large (moire lattice) as well as at the atomic scale 
[Fig. I^b)]. If a non-reactive tip is used (CO-terminated Ir 
tip) [Fig.j^c)], then imaging of the structure is different; at 
large distances only ATOP places can be clearly resolved; 
at shorter distance the Pauli repulsion between graphene 
and CO leads to the appearance of the honeycomb struc¬ 
ture formed from carbon atoms; inversion of the contrast is 
observed only at the large scale of the moire unit cell and 
not at the atomic scale, which is fully consistent with the 
Af{z) curves presented in Fig.j^c). 

In the similar 3D AFM experiments on the graphene/ 
Ru(OOOl) system, the mechanical properties of the graphene 
nanodomes, which are of the height of 1.1 A were in¬ 
vestigated ||73]| . The CFS AFM images of this system 


demonstrate the vertical reversible deformation of the hill 
places of the moire structure. The systematic measure¬ 
ments and comparison with theory allow to estimate the 
stiffness of the graphene nanomembrane to the value of 
kdome = 43.6 ± 0.5N/M with the resonance frequency of 
~ 2 THz and it was suggested that such nanodomes can be 
used as nanoelectromechanical resonators. 

Graphene moire structures on the lattice-mismatched 
metal surfaces demonstrate also the local variation of the 
chemical potential. For example, on the strongly corru¬ 
gated graphene/Rh(l 11) and graphene/Ru(0001) surfaces 
the periodic modulation of the local work function is 
0.22 eV ll74l and 0.52 eV ll75l . respectively, as obtained 
in the spectroscopic measurements of the image potential 
states with STM above different places of graphene moire. 
These results were confirmed by DFT calculations for 
the strongly interacting graphene/Ru(0001) system with 
large corrugation llT^ and explained by the different lo¬ 
cal interaction strength for topographically high ATOP 
and low FCC and HCP positions. For the weakly inter¬ 
acting graphene on Ir(lll), which has small corrugation 
of « 0.3 A with the mean distance between graphene and 
metal surface of 3.3 A, the variation of the local work func¬ 
tion in graphene moire is « 100 meV as measured from 
I{z) spectroscopy data in STM ifTTl . At the same time the 
KPFM measurements give a lower value of 35 meV ll68ll 
that can be connected with the underestimation of this 
value in KPFM measurements on the nm-scale. DFT cal¬ 
culations for this system yield value of 56 meV. 

5 Properties of graphene-metal-based systems: 
adsorption and intercalation Adsorption of graphene 
on the metal surfaces modifies the electronic spectrum of 
graphene-related valence band states around Ep. There 
are several ways allowing to tailor the properties of the 
graphene/metal interface with the aim to control electronic 
and magnetic properties of a graphene layer. 

Adsorption of different species has the aim to tailor 
the doping level of graphene as well as, in the case of the 
lattice-mismatched graphene/metal interfaces, to create the 
ordered arrays of molecules I178ll79ll80ll81ll82ll83l or metal¬ 
lic clusters 084II85I[8'6II87II88II89II where properties of the 
whole system can be modelled on the basis of the single el¬ 
ement. As proposed, such systems can be used in the future 
catalysis elements or a parts of the future information stor¬ 
age devices. For example, it was shown that the existence 
of the lateral dipole moment in the graphene/Ru(0001) 
moire structure leads to the site-selective adsorption of 
FePc and pentacene molecules on this structure and for¬ 
mation of the ordered molecular arrays [Fig. (upper 
panel)] Gl. This effect was attributed to the variation 
of the local work function along the graphene moire as 
discussed earlier. 

Adsorption of metallic clusters on top of the graphene/ 
metal moire structure was studied in a series of experi¬ 
mental and theoretical works. It was shown that for the 
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Figure 6 (a) Constant height Af{x,y) maps 

(3.5 X 2.9nm^) extracted from the 3D data set for 
graphene/Ir(lll) at different relative distances between 
tip and sample (left-hand side). The atomically resolved 
Af{z) curves are shown in the right-hand panel for dif¬ 
ferent places marked in the STM image (inset), (b) and (c) 
demonstrate the difference in the AFM and AFS imaging 
of graphene/lr(l 11) with metallic and CO-terminated tip, 
respectively. Data in (b) and (c) are taken from Ref. ca 
with permission. 


strongly corrugated graphene/Ru(0001) system adsorption 
of Pt favours to the FCC places of moire [Fig. [7] (lower 
panel)] |[85l . where as in the case of the weakly buck¬ 
led graphene on Ir(lll) adsorption of Pt and Ir clusters 
appears in HCP places 18411901 . In the latter case adsorp¬ 
tion of Ir clusters leads to the local rehybridization of 




Figure 7 Upper panel; STM images of FePc molecules 
adsorbed on graphene/Ru(0001) at varying coverage. Data 
are taken from Ref. GHl with permission. Lower panel: (a) 
STM topographic image of Pt NCs on graphene/ Ru(OOOl). 
(b) Zoom in the rectangular area in (a), (c) Line profile of 
three clusters along the purple arrow in (a), showing typi¬ 
cal heights of 1.27, 0.5, and 0.88 nm. Data are taken from 
Ref. Il85]l with permission. 


carbon orbitals from sp^ to sp^ m that was confirmed 
later in the photoemission experiments, where anisotropy 
of the Dirac cone was observed around the K point of 
the graphene Brillouine zone and the increasing of the 
width of the replica-bands-induced mini-gaps compared to 
graphene/Ir(l 11) lIMl . The further modification of this sys¬ 
tem by doping with Na atoms leads to the drastic changes 
in the spectram of graphene and opening of a large gap of 
740 meV directly at the Dirac point of graphene Il88ll . 

As has been discussed earlier the adsorption of a 
graphene layer on metal leads to the strong modification of 
the spectrum of the graphene-derived states in the vicinity 
of Ep: linear dispersion is altered via renormalisation of 
the effective mass of the carriers as well as the symmetry 
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Figure 8 STM images of (a) the graphene/Ni/Ir(lll) 
{Ut = 50mV, It = 35nA) and graphene/Cu/Ir( 111) 

(Ut = 300mV, It = 1.6nA) ll94ll intercalation-like sys¬ 
tems. Data are reproduced with permission. 


violation for the carbon sublattices in graphene might open 
the energy gap at the Dirac point. The most promising way 
to breake this “strong” interaction between graphene and 
metal is intercalation of different species, like the atoms 
of metals, gases or big molecules. This method was ini¬ 
tially developed for graphite crystals where different sub¬ 
stances were placed (intercalated) between carbon layers 
that changes the properties of graphite substantially 1921 . 

Recently the effect of intercalation of Co, Ni and Cu in 
graphene/Ir(lll) was studied II93II94||951 [Fig. [^a,b)]. As 
was found, in all cases the pseudomorphic layers of inter¬ 
calated metal (Co, Ni or Cu) were formed at the interface 
between graphene and metal. 

In case of the graphene/Ni(Co)/Ir(l 11) system II93II95I 
[Fig. [^a)] the crystallographic structure of graphene was 
found to be similar to the one for the strongly buckled 
graphene layer on Ru(OOOl); the closest distance between 
graphene and Ni (Co) layer is 1.94 A (2.02 A) and corru¬ 
gation is 1.51 A (1.27 A) compared to the same values of 
2.195A and 1.195A for graphene/Ru(0001) lIMI . It was 
found that intercalation of Ni leads to the complete redistri¬ 
bution of the valence band states of graphene and graphene 
behaves as “strongly” interacting with metal and its proper¬ 
ties (doping level as well as the formation of interface hy¬ 
brid states) are determined by the FCC and HCP “strongly” 
interacting places of the graphene/Ni/Ir(l 11) system. In the 
consequent spin-polarised STM experiments with ferro¬ 
magnetic scanning tip and in the presence of the magnetic 
field it was found that the Co islands have out-of-plane 
magnetic anisotropy (similar to Co/Ir(lll) |l96l); however 
the coercive field of the graphene/Co/Ir( 111) was found to 
be higher than ~ 4.5 T, which is much larger compared 
to the one of « 3T for Co/Ir(lll) 19^ . Comparison of 
the experimental results with DFT calculations shows that 
FCC and HCP carbon places which are closer to the Co 
layer are coupled antiferromagnetically to Co with mag¬ 
netic moment of — 1.36pB/graphene-unit-cell, whereas for 
the ATOP places, where interaction between graphene and 
Co is weaker, the coupling between Co and graphene is fer¬ 


romagnetic with magnetic moment of -fO.Id/rs/graphene- 
unit-cell. 

The situation for graphene/Cu/Ir(l 11) is similar from 
the microscopy point of view ll94l [Fig.[^b)1. The pseudo¬ 
morphic growth of Cu was found in STM experiments. At 
the same time a graphene layer on Cu/Ir(l 11) becomes less 
corrugated with the height difference of 0.229 A comapred 
to 0.307 A for graphene/Ir(l 11). Also the distance between 
graphene layer and metallic substrate is reduced after Cu 
intercalation from 3.581 A to 3.122 A for the ATOP po¬ 
sitions. It is interesting to note that despite the relatively 
“weak” interaction between graphene and Cu/Ir(lll) the 
electronic structure of graphene is significantly modified. 
Firstly, with respect to STM imaging, the inversion of the 
imaging contrast, characteristic for graphene/Ir(lll), was 
not detected and graphene/Cu/Ir(lll) is always imaged in 
the true contrast; topographically highest places (ATOP) 
are imaged as bright spots that can be explained by simi¬ 
lar carbon-projected local DOS around Ep. The electronic 
spectrum of graphene has an energy gap directly at Ep 
(similar to graphene on Cu(l 11), Cu/Ni(l 11), Ag/Ni(l 11), 
Au/Ni(lll) II97II98I 1 that was explained via reduction of 
the local symmetry for two carbon atoms in the graphene 
unit cell as a result of hybridization of Cu 3d states with 
different local character of the same interface Cu atom with 
Pz orbitals of different carbon atoms in the unit cell li94l . 

Beyond 3d metals, intercalation of Eu and Cs between 
graphene and Ir(lll) has been studied recently, reveal¬ 
ing interesting details with respect to the interfacial pro¬ 
cess Il99lll00iri01l . For the initial intercalation of Eu a for¬ 
mation of compact islands and stripes quantized in size by 
the moire unit cell are observed li99l . supporting the previ¬ 
ous reports on the growth of Ni underneath graphene on 
Rh(lll) II102I . Upon intercalation of larger amounts of 
Eu the structures grow in size yielding a pattern consist¬ 
ing of large islands and stripes divided by narrow chan¬ 
nels. This growth behaviour is attributed both to the mod¬ 
ulation of bonding strength between graphene and Ir(l 11) 
within the moire unit cell and to the strain in the graphene 
layer generated during cooling due to the difference in the 
thermal expansion coefficients of graphene and Ir. At the 
atomic level Eu forms a p{2 x 2) superstructure with re¬ 
spect to graphene, being the lowest energy configuration 
for graphene/Eu/Ir(l 11) as confirmed by the DET calcu¬ 
lations. Eurther photoemission studies showed, that partial 
intercalation of Eu or Cs leads to a strong spatial modu¬ 
lation of the graphene doping level, graphene becomes n- 
doped above the intercalated patches whereas the noninter- 
calated regions remain p-doped ifTOOl . 

Metals being the most popular intercalation materials 
are complimented by semiconductors, such as Si ifirni 
IMl, or molecular species, including hydrogen QMl, 
water II106II107I . oxygen II1081110911 1101 . carbon monox¬ 
ide cm or even Ceo II112II113I . Intercalation of these 
materials mostly aims at the electronic decoupling of the 
graphene layer form the underlying metal support. Sili- 
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con intercalated between graphene and Ru(OOOl) has been 
shown to yield well-ordered 2D patches and the disap¬ 
pearance of the strongly corrugated moire structure 111 141 . 
The latter observation coupled to the results of the cor¬ 
responding photoemission studies were considered as in¬ 
dicative for the efficient decoupling of the graphene layer 
from the Ru substrate. This method has recently been fur¬ 
ther expanded by the oxidation of the interfacial metal 
silicide layer, resulting in an insulating Si 02 film sepa¬ 
rating graphene from the metal as confirmed by photoe¬ 
mission 110411 . However, the investigations of the local 
crystallographic structure and electronic structure proper¬ 
ties of this system are still missing, thus making it difficult 
to predict if this method can be used to efficiently decouple 
graphene nanostructures from the metal substrates. 

The effective decoupling of graphene from metallic 
support can be achieved via intercalation of oxygen or 
CO 0115lll08lll09lllllllll0l . These species intercalate in 
graphene/Ir(l 11) either in the form of atomic oxygen form¬ 
ing p{2 X l)-0 layer II 1081 or CO-molecules with (3v/3 x 
3v/3)i?30° structure 01 1 111 underneath graphene. Graphene 
in these systems is fully decoupled and it is p-doped with 
the position of the Dirac point at E — Ep = 0.64 eV and 
E — Ep — 0.6 eV, respectively. For graphene/Ru(0001) 
the results are contradictory 011611117lll 101 . These works 
show that oxygen can be placed underneath graphene as 
conformed by STM and thermal-programmed-desorption 
experiments and forms (2 x 2) or (2 x 1) superstructure be¬ 
tween graphene and Ru(OOOl). However, STS experiments 
show that graphene in this system is n-doped with the posi¬ 
tion of the Dirac point at iJ—E'i;’ = —0.48 eV 01101 . that is 
opposite to the ARPES results where graphene was found 
highly p-doped with a Dirac point by « 0.5 eV higher 

EpMU- 

The interesting result was demonstrated in Ref. Una 
where intercalation of Ceo molecules in graphene/Ni(lll) 
was assumed on the basis of spectroscopic data (photoelec¬ 
tron spectroscopy and electron energy loss spectroscopy). 
Later these results were examined with STM Em¬ 
it was found that annealing of the thick layer of Ceo 
molecules leads to the desorption of most molecules, but 
also to its intercalation underneath graphene. The bias- 
dependent STM allows to discriminate between places of 
the graphene/C6o/Ni(lll) structure where intercalation 
took place and not. These measurements unambiguously 
show that Ceo intercalates via certain interfacial channels, 
which can be induced by strain relaxation at defect sites of 
graphene/Ni(lll). 

6 Graphene nanoribbons: synthesis and SPM 
studies A graphene nanoribbon (GNR) is a narrow strip 
of graphene, which morphology is determined by edge 
morphology and the width. In analogy to carbon nan¬ 
otubes, the morphology of GNRs is characterized by the 
chiral index (n, m) defining the edge translation vector 


(a) (b) 



Figure 9 (a) From top to bottom: (i) Reaction scheme from 
precursor 1 to straight N — 7 GNRs; (ii) STM image taken 
after surface-assisted C-C coupling at 200° C but before 
the final cyclodehydrogenation step, showing a polyan- 
thrylene chain (left) (Up = 1.9V, Ip = O.OSnA), and 
DFT-based simulation of the STM image (right) with par¬ 
tially overlaid model of the polymer (blue, carbon; white, 
hydrogen); (iii) Overview STM image after cyclodehy¬ 
drogenation at 400° C, showing straight N = 7 GNRs 
(RT) (Up = —3V, Ip = O.OSnA). The inset shows 
a higher-resolution STM image taken at 35 K (Up = 
— 1.5V, Ip = 0.5nA). (b) Top: Reaction scheme from 
6,1 l-dibromo-l,2,3,4-tetraphenyltriphenylene monomer 2 
to chevron-type GNRs; Bottom: Overview STM image 
of chevron-type GNRs fabricated on a Au(lll) surface 
(T = 35 K) (Up = -2V, Ip = 0.02nA). The inset shows 
a high-resolution STM image (T = 77 K) (Up = —2V, 
Ip = 0.5nA) and a DFT-based simulation of the STM im¬ 
age (greyscale) with partly overlaid molecular model of the 
ribbon (blue, carbon; white, hydrogen). Images are repro¬ 
duced from Ref. 03 with permission. 

(or, equivalently, by chiral angle, 6) and classified as either 
armchair, zigzag or chiral [see Fig.[T0|[b)]. 

The GNR electronic band structure differs from that of 
graphene due to the quantum confinement effect perpen¬ 
dicular to the GNR axis. Depending on the edge morphol¬ 
ogy and the width, GNRs may be metallic, semimetallic 
or semiconducting 020). Even from the early theoreti¬ 
cal studies employing tight binding approach it is known, 
that zigzag nanoribbons (ZGNRs) display a sharp peak in 
the electronic density of states at the Fermi level, which 
is caused by a flat band characteristic of the zigzag edge 
of graphene and leads to a net spin polarization of the 
edge ni20yi21lll22l . A small fundamental band gap opens 
up due to the antiferromagnetic coupling between the 
two spin-polarized edges. Calculations using the Hiickel 
molecular orbital method ni20y 12111122l[T2fl or a two- 
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Figure 10 (a) Representation of the gradual unzipping of one wall of a carbon nanotube to form a nanoribbon. Oxygenated 
sites are not shown. Represented from Ref. 01181 with permission, (b) A schematic drawing of an (8,1) GNR. The chiral 
vector (n, m) connecting crystallographically equivalent sites along the edge defines the edge orientation of the GNR 
(black arrow). The blue and red arrows are the projections of the (8,1) vector onto the basis vectors of the graphene lattice. 
Zigzag and armchair edges have corresponding chiral angles of 6* = 0° and 0 = 30°, respectively, whereas the (8,1) edge 
has an chiral angle of 6* = 5.8°. Lower part shows STM images of a monolayer GNR on Au(l 11) at room temperature (left) 
(Ut = 1.5V, It = lOOpA) and higher resolution STM image of a GNR at T = 7 K (right) (Ut = 200mV, It = 30pA). 
(c) Atomically-resolved STM of the terminal edge of an (8,1) GNR (Ut = 300mV, It = 60pA). (d) dl/dV spectra 
measured at different positions across the edge of GNR as shown in (c). (e) From top to bottom: edge state amplitude 
perpendicular and parallel to the edge of an (8,1) GNR and gap dependence on the width of GNRs. Data in (b-e) are 
reproduced from Ref. mu with permission. 


dimensional free massless particle Dirac’s equation with 
an effective speed of light (^ 10® m/s) 012411 125 111 261 
suggest that armchair nanoribbons (AGNRs) can be ei¬ 
ther semimetallic or semiconducting. An N-AGNR with 
N linear rows of atoms, is semimetallic if V = 3p + 2 (p 
is an integer), and semiconducting otherwise. GNRs with 
Klein edges are always semiconducting. These observa¬ 
tions are in good agreement with the data obtained with 
the first-principle approaches 0127111281 . 

The above theoretical predictions made certain GNR 
types even more attractive for applications in electronics, 
than the parent gapless material. This attraction is warmed 
up by the fact that the modern experimental techniques 
allow to create GNRs with nearly atomic precision 021. 
The problem is that the major part of theoretical predic¬ 
tions was done without consideration of an additional but 
crucial component of any working device, which is the 
metallic contacts. We have shown above, that due to the 
non-negligable degree of hybridisation between metal and 
graphene valence band states, the linear dispersion of the 
TT states in the vicinity of Ep characteristic for the free¬ 
standing graphene can be modified up to complete rear¬ 
rangement of bands. Similar effects can also be anticipated 
in GNRs. Thus, on the one hand, theoretical considera¬ 
tion of GNRs have to be revisited in order to elucidate the 


role of the substrate. On the other hand, clearly, the results 
of calculations strongly depend on the technique of calcu¬ 
lating the electronic structure and such studies have to be 
strongly linked to the experiments. 

Recently, several ways were proposed for the prepa¬ 
ration of GNRs. Among them are lithographic pattern¬ 
ing (limited by the method resolution of ~ 20 nm) mu, 
chemical sonication 11301 . gas-etching chemistry of GNRs 
produced by e-beam lithography (nanoribbons can be nar¬ 
rowed to the ^ 5nm width) 0131L cutting of a graphene 
layer along crystallographic axes by thermally activated 
metallic clusters of Fe or Ni 01321I133II134I . etc. However 
all these methods have either limited spacial resolution that 
does not allow to produce narrow nanoribbons or produce 
GNRs with edges of uncontrollable morphology that can 
lead to the uncertainty in the electronic transport properties 
of the devices built on their basis, because as was shown 
the geometry as well as defects of the GNRs edges have a 
very high impact on the electron and spin transport proper¬ 
ties of nanoribbons 01221113511136111371113801391 . 

GNRs with atomic precision can be prepared in dif¬ 
ferent ways. The first controllable synthesis was demon¬ 
strated in Ref. El, where a bottom-up approach to grow 
subnanometer wide armchair GNRs on Au(l 11) with clean 
edges was demonstrated (Fig.[^. Formation of these struc- 
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tures is performed in two thermal activation steps shown 
in panels (a) and (b) for two different organic precursors. 
In the second step the initially formed linear polymer 
chains undergo a surface-assisted cyclodehydrogenation 
and an extended fully aromatic systems are formed. The 
experimental and simulated STM images of such GNRs 
are shown in Fig. [9] Subsequent STS measurements yield 
a band gap of 2.3 ± 0.1 eV for 7-AGNR (1.4 ± 0.1 eV 
for 13-AGNRs II140I ') compared to 2.3 — 2.7eV obtained 
from quasiparticle GW calculation corrected for the im¬ 
age charge on the metallic substrate ifTSl . Parallel ARPES 
and IPES experiments give for straight 7-AGNRs and 
13-AGNRs values for the band gap of 2.8 ± 0.4 eV and 
1.6 ± 0.4 eV, respectively. Eor the chevron-type GNRs 
[Eig. |9|b)], they give a band gap of 3.1 ± 0.4eV. Also 
these measurements showed that electrons in GNRs cannot 
be considered as massless Dirac fermions as the carriers 
exhibit now a finite effective mass m* = 0.21too, where 
mo is the free electron mass. Eor the ARPES experiments, 
which require the unidirectional alignment of GNRs at 
the macroscopic scale, the similar growth procedure was 
used on the stepped Au(788) surface which consists of the 
{111} terraces of 3.83nm width II18II19I . Aligned assem¬ 
blies of both types, straight and chevron-type, GNRs can 
be grown in such case. Later, these results were examined 
within DET and many-body electron approaches CS]. 
These calculations indicate that electron transfer exists 
from any type of GNRs on Au(lll) that leads to the sur¬ 
face polarization, that is responsible for the width of the 
band gap of GNRs. The calculated energy gaps are 2.85 eV 
and 2.96 eV for 7-AGNRs and chevron-type GNR, respec¬ 
tively, that agrees well with the experimental data. 

Another approach to obtain GNRs, the so-called un¬ 
zipping of carbon nanotubes, allows flexible variation of 
GNR width, length, chirality, and substrate 01 1811 1421 
[Eig. [To} a)]. In this method the nanotubes were placed 
in an organic solution in the ultrasonic bath and then 
nanoribbons were deposited on a Au(lll) substrate and 
then their structure and electronic properties were studied 
with STM/STS Ol [Eig.[^b)]. This method allows to 
prepare GNRs of different chiralities (3.7° < 9 < 16.1°) 
and width (> 500nm) on Au(lll). As was found in STM 
the edge of GNRs is “bumped” indicating the effect of the 
GNR-metal interaction [Eig. [T0}c)]. Spectroscopic stud¬ 
ies (STS) across the edge of GNR [(c) and (d)] show the 
existence of two peaks located around Ep. The measured 
edge-state energy splitting shows a clear inverse corre¬ 
lation with GNR width [Eig. [T0|e)] and the gap values 
tend to be smaller than those observed for lithographically 
patterned GNRs (probably because of uncertainty in the 
edge structure of lithographically obtained GNRs M129II ). 
The observed spectroscopic feature (i.e. a double-peak 
in the local density of states near Ep close to the edge 
with a peak distance scaling inversely with the width of 
the GNR) was ascribed to an antiferromagnetic coupling 
between opposite edges EH, as further evidenced by 


■ > 



^ ^ 
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Figure 11 (a) Clockwise; STM overview image showing 
a free GNRs {Up = 50mV, Ip = 5pA); Eree GNR im¬ 
aged with a CO-terminated tip {Up = lOmV, Ip = 5pA); 
Constant-height high-resolution nc-AEM images of the 
middle and the zigzag end of a GNR obtained with a CO- 
terminated tip (AEM set-point offset by 30 pm); Constant- 
height nc-AEM/STM images of 22 monomer unit long 
GNR with a single missing benzene ring marked by the red 
arrow (AEM set-point offset by 48 pm). In the first image 
scale bar is 10 nm, other scale bars are 1 nm. (b) Controlled 
atomic-scale modification of the GNR reduces vibronic 
coupling; High-resolution STM image of a free GNR be¬ 
fore (Up = 50mV, Ip = 5pA) and after {Up = 50mV, 
Ip = 20pA) a bias pulse has been used to modify the left 
end of the GNR. Middle row, from left to right; Zoomed- 
in STM images of the GNR ends after the modification 
{Up = 50mV, Ip = 50pA; Up = lOmV, Ip = 2pA); 
Atomically resolved AEM image of the contacted GNR 
(AEM set-point offset by 120 pm); AEM image of the same 
ribbon as in previous image, but with the tip 80 pm closer 
to the sample. Scale bar is 0.5 nm is all images. Lower 
row; dl/dV spectra recorded at the left and right ends 
of the GNR before (blue) and after (red) the modification. 
All the images and spectra have been acquired with a CO- 
terminated tip. All data are reproduced from Ref. ED with 
permission. 

comparison with results from a Hubbard model Hamilto¬ 
nian Ml 191 [hopping integral, t = 2.7eV; on-site Coulomb 
repulsion, U = 0.5t = 1.35 eV]. These theoretical studies 
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were extended for the wider range of 0 and in the case 
of 6 = 30°, i.e. AGNR, the edge states were not present. 
The latter is in agreement with the results on AGNR/metal 
systems obtained by means of plane-wave DFT including 
semi-empirical dispersion correction 11441 . 

In agreement with the results by Tao et al. Ol, other 
experimental studies provide direct or indirect evidence for 
the presence of the edge states in supported GNRs 01451 
[US). Such edge states might be exploited for a multitude 
of spintronic applications, although it is not fully clear if 
the edge states contribute to the measured transport prop¬ 
erties of nano structures at all. To elucidate the role of 
the substrate and of edge termination, ZGNRs close to re¬ 
alistic width on (111) surfaces of Cu, Ag and Au were 
studied recently by means of DFT-D2 01471 . Although it 
was found that GNRs possess edge states, independent 
on whether they are H-free or H-terminated, they do not 
exhibit a signihcant magnetisation at the edge. These re¬ 
sults are explained by the different interaction and charge 
transfer between the GNRs and the substrates and show 
that edge magnetism in zigzag GNRs can be destroyed 
even upon deposition on a substrate which interacts weakly 
with graphene. Only in the case of H-terminated GNRs on 
Au(l 11) the interaction at the edge was found to be suffi¬ 
ciently weak to not affect the electronic and magnetic prop¬ 
erties of the edge states significantly. 

The interaction of GNRs with Au(l 11) and effect of the 
edge-termination was studied by combination of STM and 
AFM in Ref. (Ell [Fig.[g. In this work GNRs were pre¬ 
pared by the method described above 017II18I . The simul¬ 
taneous STM and AFM imaging of GNR with the tip ter¬ 
minated by CO allows to clearly identify all small changes 
in the structure of the nanoribbons. For example, the miss¬ 
ing benzene ring marked in Fig. [TTla) by the red arrow 
is clearly resolved in AFM, but not very well resolved in 
STM. Later in the same work, one of the edges of the GNR 
was modihed via removal of one of the H-atom that leads 
to the formation of the bonds between edge of GNR and 
underlying Au(lll). This place is now imaged as a dark 
(less repulsive) area at the edge of GNR [Fig. [TTlb)]. It is 
interesting to note that after such modihcation the GNRs 
cannot be moved by the STM tip (as was possible before 
modihcation), indicating the stronger interaction between 
GNR and substrate. Formation of the bond also changes the 
electronic spectrum of GNR and leads to the suppression 
of the vibronic tunnelling modes and STS spectrum shows 
now a broad peak for the left edge (H-atom is missing) of 
GNR [Fig.fTTj^b)]. Also the STS spectrum obtained on the 
intact edge (several tens nm away from the modihed edge) 
demonstrates the reduction of the intensity of the vibronic 
mode although the elastic peak remains unchanged. As was 
found in these experiments the bulk electronic structure of 
GNR is unchanged during all modihcations, but the formed 
contacts at the edges might modify the transport properties 
of GNR as the vibronic modes in the tunnelling conductiv¬ 
ity are strongly altered. 


While the interaction of nanoribbons with a noble 
metal contact is rather weak, their adsorption on metals 
like Pd and Ti, which are widely used in graphene-based 
devices, involve much stronger state mixing. The latter 
systems were considered in Ref. II148II by means of DFT. 
These calculations show that the so-called metal-induced 
gap states appear with a high density in the energy gap of 
GNRs at the metal-GNR interface contact and that they 
effectively penetrate for the large distance in GNR. These 
states can even lead to the shortening of the metal contacts 
through GNR that limits the application of the ultimately 
small GNR terminals. This work showed the importance 
of the consideration of the metal-GNR interfaces where 
intermixing of the valence band states at the interface of 
both electrodes is expected. 

7 Graphene nanoflakes and quantum dots: syn¬ 
thesis and SPM studies Graphene quantum dots 
(GQDs) are small graphene fragments (typically with the 
lateral size below 20 nm), where electronic transport is 
confined in all three spatial dimensions. GQDs can be fab¬ 
ricated by fragmentation or “cutting” of graphene sheets 
(top-down approach). Alternatively, large graphene-like 
molecules can be synthesized with well-defined molecular 
structure (bottom-up approach). 

The early work on GQDs was dominated by the inves¬ 
tigation of their physical properties. More recently, GQDs 
have been chemically modified and used for the hrst time in 
applications in the area of energy conversion, bioanalysis 
and sensors. Another exciting perspective is to use GQDs 
as spin qubits II149L The basic prerequisite is a very long 
spin-coherence time, which might exist in graphene due 
to the absence of hyperfine coupling in isotopically pure 
material and the small spin-orbit coupling. However, since 
graphene provides no natural gap, it is difficult to control 
the electron number. Moreover, the 2D sub-lattice symme¬ 
try makes the QD properties very susceptible to the atomic 
edge conhguration II14911150l unlike conventional QDs. As 
a result, chaotic Dirac billiards have been predicted 015111 
and were even claimed to be realized 01521 : i. e., the wave 
functions are assumed to be rather disordered. To achieve 
improved control of GQDs, the QD edges must be well de¬ 
fined. 

Graphene flakes or quantum dots on metals were suc¬ 
cessfully synthesised at high temperature via decomposi¬ 
tion of preadsorbed organic molecules, like hydrocarbons 
(C 2 H 4 015311231 or coronene 01541 ). or unfolding of rede¬ 
posited Cgo molecules 01550 . Fig. [T^ shows results of (a) 
STM and (b) combined AFM/KFM imaging of graphene 
flakes and QD synthesised on Ir(l 11) from C 2 H 4 by means 
of thermal programming growth (TPG) procedure 1^ . 
The next panels in the same figure show the atomically- 
resolved STM images (c) of the border between graphene 
flake and Ir(l 11 ) marked by the black frame in (a) and (d) 
single graphene QD and the corresponding edge marked 
by the white frame and the respective FFT spectrum of the 
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Figure 12 (a) and (b) 3D views of the large scale STM (160 x 160nm^, Ut = 300mV, It = O.SnA) and AFM images 
(150 X 150nm^, Af = —670mHz), respectively, of graphene flakes and QDs on Ir(lll). In (b) the topography image 
is overlaid with the results of KPFM imaging obtained simultaneously during AFM scanning, (c) Atomically resolved 
border between graphene flake and Ir(l 11) marked in (a) by black rectangle (10 x lOnm^, Ut = —50mV, It = 30nA). 
(d) STM image of the single zigzag-edge terminated GQD on Ir(lll) (20.2 x 20.2nm^, Ut = 300mV, It = 1.5nA). 
Right-hand side shows the zoomed area marked by the white rectangle (4.1 x 4.1nm^, Ut = 300mV, It = 1.5nA) and 
the corresponding FFT spectrum of the STM image of GND/Ir(l 11). 


STM image of QD where the spots corresponding to the 
reciprocal lattices of Ir and graphene are marked. As can be 
deduced from panels (a,c,d), the graphene flakes and QD 
have well-ordered structure terminated by the zigzag edges 
which are strongly bonded to the Ir(lll) surface. The 
strong surface inhomogeneity of the graphene-QD/Ir(lll) 
system allows to perform accurate AFM/KPFM measure¬ 
ments of Utcpd- Fig- [T^b) shows the 3D representa¬ 
tion of topography of the graphene-QD/Ir(l 11) system 
measured in AFM mode overlaid with the simultane¬ 
ously measured KPFM signal. The measured difference 
of Ulcpd is ~ 600 meV that is smaller than the differ¬ 
ence of work functions of Ir(lll) and graphene extracted 
from STS (1.1 ± 0.3eV) and LEEM (1.6eV) JBTI . 

This effect can be connected with the smearing effect of 
the macroscopic scanning tip with a sharp apex 0158II159I . 


Direct insight into electronic structure and properties of 
QD is provided by STS which maps out the LDOS of con¬ 
fined charge carriers and, at the same time, determines the 
shape of the QD atom by atom. Recently, several groups 
investigated in situ prepared GQDs with exclusive zigzag 
edges, which are supported by Ir(lll) II153II160 11^12^ . 
The presence of edge-states was not detected ll25l . The 
experimental results are supported by DET calculations. 
These observations are explained by a hybridisation of the 
GQDs Pz orbitals and the substrate valence band states 
(here: Ir 5^22 surface state lying at Ep), which gradually 
decreases in strength from the edge towards the center of 
the GQD. Flowever, these works, presenting the similar ex¬ 
perimental data on the confinement of the electronic states 
in GQDs, do not provide a clear answer to what extent the 
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graphene electronic states contribute to the tunnelling cur¬ 
rent compared to the contribution from the substrate. 


This problem was recently addressed in Ref. 
whose main results are presented in Fig. 13 In this work 
graphene flakes were initially produced on Ir(l 11) via TPG 
procedure and then they were decoupled from it via inter¬ 
calation of thick layer of Au. In this case the graphene- 
flakes/Au(l 11) system of very high quality is formed 
demonstrating herringbone reconstruction of Au(lll) as 
well as a moire structure on the graphene flakes due to the 
lattice mismatch of graphene and Au(lll). Compared to 
graphene/Ir(lll) where graphene flakes are bonded to the 
Ir substrate at the edges II161I . they are very mobile on the 
Au(lll) surface. This is conformed by the possibility to 
move graphene flakes across the surface by the STM tip 
with only moderate tunnelling currents (typical parameters 
are Ut = 1V and = 1 nA). 


The weak interaction between the graphene flakes 
and Au(lll) allows to carefully follow the dispersion 
of the electronic states obtained from STS experiments 
(see Fig. ES- Panel (a) shows an STM image of a single 
graphene flake on Au(l 11) and panel (b) presents the cor¬ 
responding dl/dV map measured at Ut = —75 mV of 
the area from (a) where simultaneously the Au(lll) sur¬ 
face and the graphene-flake/Au(l 11) interface are probed. 
One can clearly see that on both sides of the border the 
electronic standing waves have the same character that 
indicates the same origin of the electronic states - surface 
sate of Au(l 11). Similar spectroscopy maps were acquired 
at different bias voltages. FFT analysis of these dl/dV 
maps demonstrates the existence of two circles in the re¬ 
ciprocal space with the radii of Qgr/Au qau, which can 
be assigned to the surface state of Au(lll) on graphene- 
covered and graphene-free surfaces. The resulting FFT 
data allow to plot the respective dispersions of the elec¬ 
tronic states E{k), where qAu,gr/Au = [Fig. 13 f)]. 
The extracted effective mass is similar for both electronic 
states, m* = 0.26me. 


The similar dl/dV measurements with atomic res¬ 
olution were performed on the same graphene flake 
[Fig. [T3]^d)]. The extracted FFT pictures show additional 
features: (i) six spots corresponding to the reciprocal lattice 
of graphene, which are superposed by the spots originating 
from the reconstruction of Au(lll) and the moire struc¬ 
ture (large hexagon); (ii) ring-like features which build a 
hexagon corresponding to the (v/3 x v/3)7?30° structure 
in the real space and related to the intervalley scattering 
(qinter — UK — 2k). The obtained FFT data at different 
bias voltages allow to plot the dispersion of these elec¬ 
tronic states [red squares in Fig. [T3|f); they are plotted 
with respect to the K point of the graphene Brillouine 
zone]. The linear fit of these data gives the Fermi veloc¬ 
ity of vp = (1.1 ± 0.2) • 10® m/s and the position of the 
Dirac point at E — Ep = 0.24 ± 0.05 eV, i. e. graphene 
is p-doped as found in the earlier photoemission studies of 
the graphene/Au(l 11) system lfT62l . These measurements 


performed on the same graphene flake show the spectro¬ 
scopic features from the surface state of Au(l 11) and from 
the intervalley scattering of the graphene Dirac fermions 
allow to give an answer about different contributions in the 
imaging of graphene on metals Ill53lll60ll23ll24ll26lll63l 
[1641 . 

The similar experiments were performed on the 
GQDs/oxygen/lr(l 11) system nsi. A linear fit of the data 
gives the Fermi velocity of uf = (0.96±0.07)-10® m/s and 
the position of the Dirac point at E — Ep = 0.64±0.07 eV, 
which agrees well with the data obtained in ARPES exper¬ 
iments in the same and similar works 1163111091 . 

Preparation of graphene flakes or GNDs on noble met¬ 
als might lead to the formation of highly strained graphene 
nonobubbles, as it was demonstrated in Ref ll65l for 
graphene-flakes/Pt(l 11). Such structures have unique elec¬ 
tronic properties as was demonstrated in STS measure¬ 
ments. The STS dl/dV spectra measured across such 
graphene nanobubble show several spectral features, which 
can be assigned to the Landau levels in graphene. The na¬ 
ture of this effect is the existence of the so-called pseudo 
magnetic field (of more than SOOT) originating from the 
strain appearing in the graphene nanobubbles during prepa¬ 
ration (different expansion coefficients of graphene and 
Pt). The similar STS signatures were found for graphene 
nanobubbles formed by the intercalation of oxygen in 
graphene/Ru(0001) 11661 . finally giving a chance to form 
the regular arrays of the strain-engineered graphene-based 
structures. 

8 Conclusions The present review demonstrates the 
powerful potential of the modern SPM methods in con¬ 
junction with the corresponding spectroscopy for the in¬ 
vestigation of the graphene-based systems on metallic sur¬ 
faces. These recent experimental works are devoted to the 
studies of the graphene-metal interfaces and nanostructures 
where different phenomena were observed: localisation of 
the electronic states in graphene nanoribbons and GND, 
edge electron scattering in different geometries, demon¬ 
stration of the Landau levels in strained graphene nanos¬ 
tructures, etc. For simple graphene-flakes on metals it was 
shown that it is possible to discriminate between graphene- 
and metal-related contributions in the electronic structure 
of the system at the nanoscale. All these results are found 
to be in very good agreement with state-of-the-art DFT cal¬ 
culations that allow to model more complicated graphene- 
based nanostructures. Such an approach permits to study 
the more realistic graphene-metal systems which will find 
a broad utilisation in future electron- and spin-transport ap¬ 
plications and devices. 
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Figure 13 (a) Topography of an approximately 400 x 160 nm^ large flake used for dl/dV mappings, (b) dl/dV map 
across the rim of a large graphene flake showing the quasiparticle interferences of the Au surface state {Ut = —75mV, 
It = InA). (c) Fast Fourier transform of a 109 x 109 nm^ large area partially depicted in (b) at selected bias voltages, 
(d) dl/dV map on the graphene flake with atomic resolution (Ut = —20mV, It = InA). (e) Fast Fourier transform of a 
54 X 54 nm^ map partially depicted in (d) showing a rich structure including atomic, moire, and herringbone features as well 
as intervalley scattering features. The right column shows magnifications of the intervalley scattering rings (3.5 x 3.5 nm“^) 
at different bias voltages, (f) The Au(lll) surface state dispersion for pristine Au(lll) (cyan diamonds) as well as for 
graphene/Au( 111) (blue triangles). Dispersion relation of the graphene electrons determined from the intervalley scattering 
(red squares) and corresponding fit (red line) including uncertainty (red dotted lines) determined from a series of atomically 
resolved constant-energy mappings. The k values are plotted with respect to the T-point in the case of the surface state 
and with respect to the AT-point in case of graphene. Data are reproduced from Ref. Il26l with permission. 
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